Critical behavior of repulsive linear fc-mers on triangular lattices 
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Monte Carlo (MC) simulations and finite-size scaling analysis have been carried out to study the critical 
behavior in a submonolayer two-dimensional gas of repulsive linear fc-mers on a triangular lattice at coverage 
k/(2k+ 1). A low-temperature ordered phase, characterized by a repetition of alternating files of adsorbed k- 
mers separated by k + I adjacent empty sites, is separated from the disordered state by a order-disorder phase 
transition occurring at a finite critical temperature, T^. The MC technique was combined with the recently re- 
ported Free Energy Minimization Criterion Approach (FEMCA), [F. Roma et al., Phys. Rev. B, 68, 205407, 
(2003)], to predict the dependence of the critical temperature of the order-disorder transformation. The de- 
pendence on k of the transition temperature, Tc{k), observed in MC is in qualitative agreement with FEMCA. 
In addition, an accurate determination of the critical exponents has been obtained for adsorbate sizes ranging 
between k = I and k = 3. For k> I, the results reveal that the system does not belong to the universality class 
of the two-dimensional Potts model with q = 3 (k= 1, monomers). Based on symmetry concepts, we suggested 
that the behavior observed for A: = 1 , 2 and 3 could be generalized to include larger particle sizes (k>2). 
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I. INTRODUCTION 



Two previous articles^'", referred to as Papers I and II, re- 
spectively, were devoted to the study of repulsively interact- 
ing A:-mers on square lattices at half coverage. In Paper I, 
Monte Carlo (MC) simulations were used to deal with such 
problem. A {2k x 2) ordered phase, characterized by a rep- 
etition of alternating files of adsorbed A:-mers separated by k 
adjacent empty sites, was found. This ordered phase is sepa- 
rated from the disordered state by a order-disorder phase tran- 
sition occurring at a critical temperature Tc. The depen- 
dence on the particle size shown an intriguing behavior, with 
a pronounced minimum at A: = 2. In addition, two theoreti- 
cal approaches. Detailed Mean-Field Approximation (DMFA) 
and Free Energy Minimization Criterion Approach (FEMCA), 
were combined with Monte Carlo simulations to obtain the 
critical temperature of the order-disorder transformation. Pre- 
dictions from FEMCA, based on a free energy minimization 
criterion, shown a remarkable qualitative agreement with the 
simulation data and allowed us to interpret the physical mean- 
ing of the mechanisms underlying the observed transitions. 
Paper II was a step further, analyzing the universality class 
of the phase transition. For this purpose, an extensive work 
of MC simulations and finite-size scaling analysis was car- 
ried out. Based on the strong axial anisotropy of the low- 
temperature phase for k>2, a new order parameter measur- 
ing the orientation of the admolecules was introduced. Tak- 
ing advantage of its definition, an accurate determination of 
the critical exponents was obtained for adsorbate sizes rang- 
ing between k = 2 and k — A. The evaluation of the complete 
set of static critical exponents, a, J and V for different 
molecular sizes, shown that: for linear A;-mers with k>2, this 
phase transition does not belong to the universality class of 
the two-dimensional Ising model {k = 1). The main source 
for this behavior was associated to the breaking of the orien- 
tational symmetry occurring for k>2, which does not occur 



for k = 1. Moreover, the critical exponents reported in Paper 
II do not correspond to a known universality class, accord- 
ing to the current classification of order-disorder transitions 
on surfaces given by M. Schick^. The aim of this paper is the 
study of repulsive straight A;-mers on triangular lattices using 
the same techniques developed in Papers I and II. 

The behavior of interacting dimers has also been analyzed 
by using graph theoryi. In Ref4., a Sierpinski gasket with a 
fractal dimension of In3/ln2 was used as substrate. How- 
ever, the technique can be extend to other fractal graphs with 
low "ramification degree". With respect to triangular lat- 
tices, leading contributions have been published by Hock and 
McQuistan^ and Phares and Wunderlich^. The thermodynam- 
ics of non-interacting dimers was obtained from the knowl- 
edge of the occupational degeneracy of dimers partially cov- 
ering the lattice. Later, Rzysko and Borowkc^ studied the in- 
teresting problem of heteronuclear dimers consisting of dif- 
ferent segments A and B adsorbed on square and triangular 
lattices. The authors considered models with attractive B B 
and A — B nearest-neighbor energies and variable A — A en- 
ergy (attractive as well as repulsive). A rich variety of phase 
transitions was reported along with a detailed discussion about 
critical exponents and universality class. In the case of a tri- 
angular lattice all systems belong to the class of universality 
of the two-dimensional Potts model with q — 3. 

In previous work we have studied the adsorption thermo- 
dynamics of attractive and repulsive dimers on honeycomb 
and triangular lattices^. From the study of adsorption proper- 
ties such as adsorption isotherm, heat of adsorption and con- 
figurational entropy, two different ordered phases occurring 
in the adlayer were reported for repulsive couplings and low 
temperatures. Namely, (1) the low coverage ordered phase 
(LCOP), with 2/5 (5/9) of the sites occupied for a triangular 
(honeycomb) lattice and (2) the high coverage ordered phase 
(HCOP), with 2/3 of the sites filled for both lattices. In the 
particular case of triangular lattices, the LCOP resembles the 
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(4 X 2) ordered structure of repulsive dimers on square lat- 
tices at half coverage. As in the last case, from the LCOP 
appearing in dimers at = 2/5, we can predict the existence 
of a ordered structure for repulsive straight fc-mers on trian- 
gular lattices at coverage kj {2k + 1). A snapshot of this low- 
temperature phase, characterized by a repetition of alternating 
files of adsorbed ^-mers separated by A:+ 1 adjacent empty 
sites, is shown in Fig. 1 for k — 2. 

In this context, the scope of the present work is: (/) to calcu- 
late, via MC simulationsiiSiii, finite-size scaling analysis!^ 
and FEMCA, the critical temperature as a function of the 
size of the adsorbed molecules for repulsive linear A:-mers ad- 
sorbed on triangular lattices at = k/{2k + 1 ); and (//) to de- 
termine critical exponents and universality. The outline of the 
paper is as follows: In Section 2 we describe the lattice-gas 
model and the simulation scheme. In Section 3 we present the 
results. Finally, the general conclusions are given in Section 
4. 



II. LATTICE-GAS MODEL AND MONTE CARLO 
SIMULATION SCHEME 

A. The model 

In this section we describe the lattice-gas model for the ad- 
sorption of particles with multisite occupancy in the mono- 
layer regime. We consider the adsorption of homonuclear k- 
mers on triangular lattices. The adsorbate molecules are as- 
sumed to be composed by k identical units in a linear array 
with constant bond length equal to the lattice constant a. The 
k-meis can only adsorb flat on the surface occupying k lat- 
tice sites. The surface is represented as a triangular array of 
M = LxL adsorptive sites, where L is the linear size of the ar- 
ray. In order to describe a system of ^-mers adsorbed on M 
sites at a given temperature T, let us introduce the occupation 
variable c, which can take the values c,- = or 1, if the site / 
is empty or occupied by a A:-mer unit, respectively. The A:-mer 
retain its structure upon adsorption, desorption and diffusion. 
The Hamiltonian of the system is given by. 



-N{k^ l)w + e„£ci 



(1) 



where w is the nearest-neighbor (NN) interaction constant 
which is assumed to be repulsive (positive), (/, /') represents 
pairs of NN sites and £„ is the energy of adsorption of one 
given surface site. The term N{k— l)w is subtracted in Eq. 
Q since the summation over all the pairs of NN sites overesti- 
mates the total energy by including N{k — I) bonds belonging 
to the adsorbed k-mers. Finally, £„ is set equal zero, without 
any lost of generality. 



B. Order parameter and Monte Carlo method 

In order to study the order-disorder phase transition occur- 
ring in the adsorbate, it is convenient to define a related order 



parameter. For k = 1, a standard geometrical order parame- 
ter can be built^'*. For this purpose, the lattice is separated in 
three intersecting sublattices, p = 1,2,3, and for each one of 
them the functions: 
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are obtained. Then the order parameter (p of the system is 
given by: 
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where 
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and the 5 takes the value 1 for a cycUc permutation of 
subindexes (1,2,3) and the value otherwise. 

As in Paper II, an order parameter measuring the orientation 
of the admolecules in the ordered structure will be defined 
for k > I. Fig. 1 shows a snapshot corresponding to one of 
the possible configurations of the low-temperature phase for 
dimers at = 2/5. As it can be observed, the structure adopts 
an orientation along one of the lattice's axis. Taking advantage 
of this property, we define dk as: 
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where Nx{x= 1,2,3) represents the number of A:-mers aligned 
along one of the three axis of the lattice and = A^i +N2+N3. 
It is worth to notice that, due to the low-temperature phase is 
isotropic for monomers, 5k is restricted to > 2. 

When the system is disordered {T > T^), the three orien- 
tations are equivalents and 5ii is zero. As the temperature is 
decreased below Tc, the A:-mers align along one direction and 
5jt is different from zero (being 5^ = I at T = 0). Then, 5^ ap- 
pears as a good order parameter, evidencing the order-disorder 
phase transition. 

As it is well-known for k = I, the system belongs to the 
universality class of the three-state Potts mode l?i^?i'^i^? , being 
V = 5/6, j3 = 1/9, 7 = 13/9 and a 1 /3. In this case, the 
phase transition is accomplished by a breaking of the trans- 
lational symmetry (the coverage of the sublattices is different 
for T < Tc). However, an additional breaking of the orien- 
tational symmetry occurs for k>2. Consequently, a change 
in the universality class is expected for linear molecules with 
k>2, with respect to the case of monomers. 

In order to study the critical behavior of the system, we 
have used an efficient exchange MC method'^'^^ and finite- 
size scaling analysis2iiSi22. As in Ref«i^ we build a com- 
pound system which consists of m non-interacting replicas of 
the system concerned. The m-th replica is associated with the 
temperature r„, [or j3,„ = l/{kBTm), being kg the Boltzmann 
constant]. In other words, each replica is in contact with its 
own heat bath having different temperature. Under these con- 
ditions, the algorithm to carry out the simulation process is the 
following: 
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1) The compound system of m replicas is generated. For 
this purpose, each repHca is simulated simultaneously 
and independently as canonical ensemble for ni MC 
steps by using a standard importance sampling MC 
method^ '"-^'. In order to determine the set of tempera- 
tures, {j3m}, we set the highest temperature, T,„ax{^max), 
in the high temperature phase where relaxation (corre- 
lation) time is expected to be very short and there ex- 
ists only one minimum in the free energy space. On the 
other hand, the lowest temperature, r,„,„ (j3,„,„ ) , is some- 
where in the low-temperature phase whose properties 
we are interested in. Finally, the difference between two 
consecutive temperatures is set as {T,ru,x — Tmm) / ('w — 
1) (equally spaced temperatures). 

2) Interchange vacancy-particle and diffusional 
relaxations^ The procedure is as follows: 

2.1) One of the m replicas is randomly selected; 

2.2) a A:-mer and a linear A:-uple of empty sites, both 
belonging to the replica chosen in 2.1), are ran- 
domly selected and their positions are established. 
Then, an attempt is made to interchange its occu- 
pancy state with probability given by the Metropo- 
lis rule^': 

P = min{l,exp(-A„A//)} (6) 

where AH = Hf — Hi is the difference between the 
Hamiltonians of the final and initial states; 

2.3) a fc-mer is randomly selected. Then, a dis- 
placement to nearest neighbor positions is at- 
tempted (following the Metropolis scheme), by ei- 
ther jumps along the A;-mer axis or reptation by ro- 
tation around a unity of the A:-mer This procedure 
(diffusional relaxation) must be allowed in order 
to reach equilibrium in a reasonable time. 

3) Exchange of two configurations X„, and X,,/, cor- 
responding to the m-th and w'-th replicas, respec- 
tively, is tried and accepted with the probability 
W (X„,,j3,„|X„,/,/3„,'). In general, the probability of ex- 
changing configurations of the m-th and m'-th replicas 
is given b}*™, 

W(X,„,A„|X,„,,^„,0 = {e,p(_A) trltl 

where A = (/3,„ - A,,') [i/(X/„) -i/(X,„)]. As in Ref A 
we restrict the replica-exchange to the case m ^ m + 1 . 

4) Repeat from step 2) m x M times. This is the elemen- 
tary step in the simulation process or Monte Carlo step 
(MCS). 

The procedure l)-4) is repeated for all lattice's sizes. For 
each size, the equilibrium state can be well reproduced after 
discarding the first n2 MCS. Then, averages are taken over 



riMcs successive MCS. The canonical expectation value of a 
physical quantity A is obtained in the usual way as follows: 

1 nMCS 

{A)p^ = l^A[X„{t)] (8) 

All calculations were carried out using the parallel cluster 
BACO of Universidad Nacional de San Luis, Argentina. This 
facility consists of 60 PCs each with a 3.0 MHz Pentium-4 
processors. The values of the parameters used in the simulated 
tempering runs are shown in Table I. 

The internal energy per lattice site, u, is obtained as simple 
averages: 

u^^{H)^. (9) 
The specific heat C is sampled from energy fluctuations. 

The quantities related with the order parameter, such as the 
susceptibility ™d the reduced fourth-order cumulant Ul 
introduced by Binde»^, can be calculated as: 

X = ^[{z^)t-{z)\] (11) 



t/.(r) = i-|| (12) 

where z(= (p or 5k) represents the order parameter and the 
thermal average (...)7', in all the quantities, means the time 
average throughout the MC simulation. 



III. RESULTS 

The calculations were developed for linear A:-mers, with k 
ranging between 1 and 3, on a triangular L x L lattice at cov- 
erage k/ {2k+ 1). In addition, conventional periodic boundary 
conditions were considered. Note, however, that the choice of 
appropriate linear dimensions L has to be done in such away 
that the ordered structures are not disturbed. Thus, for k=\, 
the effect of finite size was investigated by examining lattices 
varying from L = 30 to L = 120. In the case of A; = 2 and 
/t==3, L= 10, 15,20,25,30 and L= 14,21,28,35 were used, 
respectively, with an effort reaching almost the limits of our 
computational capabilities. 

The critical temperatures have been estimated from the 
plots of the reduced four-order cumulants Ul{T) plotted ver- 
sus kgT /w for several lattice sizes. In the vicinity of the crit- 
ical point, cumulants show a strong dependence on the sys- 
tem size. However, at the critical point the cumulants adopt 
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a nontrivial value U*\ irrespective of system sizes in the scal- 
ing limit. Thus, plotting Ul{T) for different linear dimensions 
yields an intersection point U*, which gives an accurate esti- 
mation of the critical temperature in the infinite system. 

Hereafter we discuss the behavior of the critical tempera- 
ture as a function of the size k, for adsorbed A:-mers at mono- 
layer Fig. 2 illustrates the reduced four-order cumulants 
Ul{T) plotted versus kgT /w for several lattice sizes and k~l. 
From their intersections one gets the estimation of the critical 
temperature. The value obtained for the critical temperature 
was fcB7;./w = 0.3354(1), being f/* =0.617(1). Our calcula- 
tion of the critical temperature is in well agreement with the 
transfer-matrix result of Kinzel and Schick^^. With respect to 
Monte Carlo calculations, Metcalf^^ in the 70's and Chin and 
Landau^^ in the 80's obtained similar values of the critical 
temperature. Namely, kgTc/w w 0.35. Due to the lattice sizes 
and the number of MCS used in this contribution, our esti- 
mate of Tc is expected to be more accurate than those reported 
previously. 

Figs. 3 and 4 show the results of Ul{T) for k — 2 and k — 3, 
respectively. In the case k — 2, the value estimated for kgT^./w 
[= 0.2338(1)] agrees very well with previous determinations 
reported in the literature^. In Ref.^, a value ksTc/w ~ 0.23 
was obtained from the inflection on the function s{T), being 
s{T) the configurational entropy of the adlayer as a function 
of the temperature. With respect to A: = 3, the value obtained 
for ksTc/w [= 0.341(1)] is reported for the first time in the 
literature. 

The values of the critical temperature, which are collected 
in Table II for k — 1 — 3, present a non-trivial behavior as a 
function of the particle size k. To understand the dependence 
of Tc{k) on k, it is convenient to appeal to the recently reported 
FEMCA^ In this theoretical framework, Tc depends on the 
ratio of the energy per lattice site and entropy per lattice site 
differences, u„ — mq and s„ — sq, respectively, between a fully 
disordered state {T oo) and the ground state {T 0), 
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In this case, the mean energy and the entropy for the ordered 
state are uq — sq— 0. Then, the critical temperature depends 
on the mean energy and the entropy of the disordered state. Moo 
can be calculated from the quasi-chemical approximation^^ 
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being c the lattice connectivity. In addition, 

A = {c-2)k + 2, 



a = 



0(1-0) 



2k 



' k-\ ^ 
y k , 



+ ^7 



(14) 



(15) 



(16) 



and 



k-1 



— [l-exp(-w/feBr)] 0(1-0) j 

(17) 



With respect to the entropy. 



and 
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Soo//tB = ln3-(2/3)ln2 forA:=l, 
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for k>2. (19) 



In the last equation, the entropy was obtained from the ex- 
tension to higher dimensions of the exact configurational fac- 
tor of linear chains adsorbed on one-dimensional lattices ». 
Under these considerations, the critical temperature results, 

ksTcik) Moo/w 
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for ;t = 1 
for A: > 2 



Fig. 5 presents the theoretical predictions obtained from 
FEMCA (Eq. I20> for the critical temperature as a function of 
the size k. For high values of k, kBTc{k)/w increases mono- 
tonically as k is increased (see inset). In the other extreme, the 
theoretical curve shows a marked break at = 2, which is in 
good qualitative agreement with the minimum obtained from 
Monte Carlo simulations. 

The theoretical approach allows us to interpret the physical 
meaning of the main features of the critical temperature. In 
fact, Eq. il3\ shows that ksTc/w depends on the mean energy 
and the entropy of the disordered state. In first term, we will 
analyze the entropy. When one pass from {k = 1,0 = 1/3) 
to (^ = 2,0 = 2/5), the number of entities per lattice site, 
«£, and for this reason, the number of accessible states, di- 
minishes [in general, ne{k) oc 1 /{2k+ 1)]. On the other hand, 
a new degree of freedom appears for the adsorbed particles 
whose size is ^ > 2: dimers, trimers, etc. can rotate on the 
lattice. The new accessible states for dimers due to possible 
rotations compensate the diminution in the number of enti- 
ties with respect to monomers. Consequently, the variation 
in the entropy between monomers and dimers is small. Thus, 
s„{k = 1,0 = 1/3)= 0.6365 [see Eq. (d}] and s^{k = 2,0 = 
2/5) « 0.6339*^ for monomers and dimers adsorbed in trian- 
i/2gular lattices, respectively. However, due to the new degree 
of freedom appears for k>2, the diminution in the entropy is 
reestablished (associated to the diminution in the number of 
entities) and Soo diminishes as k is increased for k> 2. 
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With respect to the mean energy, the approximated so- 
lution given in Eq. J14l i shows that Uoo{k) increases until 
an asymptotic value for higher fe's: Moo(1) = 0.333, Moo(2) = 
0.357, Moo(3) =0.368, M„o(4) =0.375, •••,M„oH =0.40. 

Now we can interpret the two different regimes observed 
in the curve of Fig. 5. While Uoo{k) changes slightly in all 
range of k, Soo{k) presents an abrupt change at A: = 2. Physi- 
cally, this change in the entropy, which is responsible for the 
marked break appearing in the curve of ksT^lk) /w versus k, is 
associated with the possibility of orientation of the adsorbed 
molecules for k>2. 

As it was shown in Figs. 3 and 4, the intersections of the 
cumulants for k — 2 and k — 3 suggest that the system suffers 
a phase transition at coverage k/ {2k +1). In order to corrob- 
orate this finding, the critical exponents will be calculated in 
the following. 

In first term, we study exhaustively the well-known lattice- 
gas of repulsive monomers. Once we know 7^ (see Fig. 2), 
the critical exponent v can be calculated from the full data 
collapsing of Ul- The results are shown in Fig. 6 (being 
t = T/Tc~l), where an excellent fit was obtained for v = 5/6. 
Given ksTc/w = 0.3354(1 ) and v = 5 /6, a, /3 and y were ob- 
tained from the collapse of the curves of C, (p and as it is 
shown in Figs. 7-9, respectively. The data scaled extremely 
well using a = 1/3, j3 = 1/9 and 7= 13/9, which corrob- 
orates that the case of A: = 1 corresponds to the universality 
class of the two-dimensional Potts model with q — 3. The val- 
ues of the critical exponents along with the intersection point 
of the cumulants are collected in Table II. 

The finite-size scaling study was carried out for A: = 2 (Figs. 
10-13) and k = 3 (data do not shown here). The resulting 
values of the critical exponents, which are listed in Table II, 
point out clearly the existence of a different universality class 
from that of the three-state Potts model. 

As can be demonstrated, the set of critical exponents for 
dimers and trimers fulfills the well-known inequalities of 
Rushbrooke^^, a + 2j8 + 7 > 2, and Josephson^**, dv + a>2 
(being d the dimension of the space), leading to independent 
controls and consistency checks of the values of all the critical 
exponents. 

The identical results (within numerical errors) obtained for 
the critical exponents corresponding to dimers and trimers 
corroborate our hypothesis that the breaking of the orienta- 
tional symmetry occurring for k>2: i) affects the nature of 
the phase transition and //) is the main source of the change 
in the universality class, with respect to the well-known two- 
dimensional three-state Potts class of universality correspond- 
ing to monomers. 



Surprisingly, the exponents obtained for repulsive A;-mers 
(fc > 1 ) on a triangular lattice at9 =k/ {2k + 1 ) are very similar 
to those characterizing the critical behavior of A^-mers (A^ > 1 ) 
on square lattices at = 1 /2^. However, we think that these 
two phase transitions belong to different universality classes. 
Some evidences support our arguments: 1) the value obtained 
for the fixed-point cumulant U* varies from square to triangu- 
lar lattices^^; and 2) the low-temperature phase have two and 
three possible orientations for square and triangular lattices, 
respectively. Then, the order parameters characterizing both 
phase transitions have different symmetries. 



IV. CONCLUSIONS 

In the present work, we have addressed the critical proper- 
ties of repulsive linear A:-mers on two-dimensional triangular 
lattices at coverage k/{2k+ 1), and shown the dependence of 
the critical temperature on the size k. The results were ob- 
tained by using MC simulations and theoretical calculations 
from FEMCA. 

The critical temperature dependence on the particle size of 
the low coverage ordered phase [9 — k/{2k+ 1)] of repul- 
sive straight particles have been reported for the first time, 
and found that dimers present the minimum value. FEMCA 
shown that an analysis of the delicate balance between the size 
dependence of entropy and energy per site allows us to inter- 
pret qualitatively the behavior of Tc versus k. 

On the other hand, the evaluation of the complete set of 
static critical exponents, a, P, y and v for different molecu- 
lar sizes, shows that for linear A^-mers with k>2, this phase 
transition does not belong to the universality class of the two- 
dimensional Potts model with q = 3. The main source for this 
behavior is the breaking of the orientational symmetry occur- 
ring for k>2, which does not occur for k= I. Moreover, the 
critical exponents reported in the present paper for A: > 2 do 
not correspond to a known universality class, according to the 
current classification of order-disorder transitions on surfaces 
given by Schicb^ 

A detailed scaling analysis shows that usual hyperscaling 
relations (inequalities of Rushbrooke and Josephson) are ful- 
filled and leads to independent controls of the values of all the 
critical exponents. 

Future efforts will be directed to obtain the critical behavior 
for other existing ordered phases in the whole range of cover- 
age. 
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V. TABLES AND FIGURE CAPTIONS 

Table I Parameters of the simulated tempering runs. 

Table II Values of Tc, a, j3, 7, V and U* (as indicated in the 
text) for k ranging from 1 to 3. In the case k = I, the 
critical exponents correspond to the well-known two- 
dimensional Potts model with q~3. 
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Snapshot of the ordered phase for dimers at 9 — 2/5. 

Curves of Ul{T) versus kgT /w, for k—\ and different 
lattice sizes as indicated. From their intersections one 
obtained ksTc/w. 

Same as Fig. 2 for dimers (k = 2). 

Same as Fig. 2 for trimers (k = 3). 

Comparison between simulated and theoretical results 
for ksTc/w vs. k. The symbology is indicated in the 
figure. 

Data collapsing for the cumulants in Fig. 2. 

a) Size dependence of the specific heat, C (with k ~ 
1) as a function of temperature, b) Data collapsing 
{CL-"!" vs. fL^/'^) for the curves in a). 

a) Size dependence of the order parameter, <p (with 
A; = 1) as a function of temperature, b) Data collaps- 
ing {(pL^I^ vs. |f |L'/^) for the curves in a). 

a) Size dependence of the susceptibility, % (with k = 
1) as a function of temperature, b) Data collapsing 
(jL^^/*' vs. tL}l^) for the curves in a). 

Data collapsing for the cumulants in Fig. 3. 

Same as Fig. 7 for dimers {k = 2). 

Same as Fig. 8 for dimers {k~2). 



Figure 13: Same as Fig. 9 for dimers (k — 2). 
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TABLET 



k 


L 


m 


«i 


«2 




Tmin 


Tmax 




30 


40 


10^ 


10^ 


5 X 10^ 


0.279 


0.392 


1 


60 


40 


10^ 


10^ 


5 X 10^ 


0.311 


0.360 




90 


40 


10^ 


10^ 


5 X 10^ 


0.290 


0.381 




120 


40 


10^ 


10^ 


5 X 10^ 


0.325 


0.346 




10 


41 


5 X 10"* 


5 X 10"* 


10^ 


0.22 


0.30 




1 


41 


5x10^ 


5x10''^ 


10^ 


0.22 


0.27 


2 


20 


61 


5 X lO'* 


5 X 10'^ 


10^ 


0.22 


0.26 




25 


101 


5x10^ 


5x10^ 




0.23 


0.25 




30 


131 


2.5 X 10-^ 


2.5 X 10^ 


5 X 10 -"^ 


0.23 


0.25 




14 


101 


1.5 X 10^ 


1.5 X 10^ 


3 X 10^ 


0.300 


0.400 




21 


61 


2.5 X 10^ 


2.5 X 10^ 


5 X 10^ 


0.320 


0.380 


3 


28 


51 


2.5 X 10^ 


2.5 X 10^ 


5 X 10^ 


0.335 


0.360 




35 


51 


10^ 


10^ 


2x 10^ 


0.337 


0.347 



TABLE II 



k 


ksTc/w 


a 




r 


V 


U* 


1 


0.3354(1) 


1/3 


1/9 


13/9 


5/6 


0.617(1) 


2 


0.2338(1) 


0.95(3) 


« 0.0005 


1.17(3) 


0.50(1) 


0.656(1) 


3 


0.341(1) 


1.00(4) 


« 0.0000 


1.20(4) 


0.50(1) 


0.656(2) 



k = 2; 9=2/5 



Fig. 1 : P. M. Pasinetti et al. 
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